InSAR bias and uncertainty due to the systematic and stochastic tropospheric delay

We quantify the bias and uncertainty of interferometric synthetic aperture radar (InSAR) displacement time series and their derivatives, the displacement velocities, by analyzing the systematic and stochastic components of the temporal variation of the tropospheric delay. The biases due to the systematic seasonal delay depend on the SAR acquisition times, whereas the uncertainties depend on the standard deviation of the random delay, the number of acquisitions, the total time span covered, and the covariance of the time series of the stochastic delay between a pixel and the reference. We study the contribution of the wet delay to the InSAR observations along the western India plate boundary using (i) Moderate Resolution Imaging Spectroradiometer precipitable water vapor, (ii) stratified tropospheric delay estimated from the ERA-I global atmospheric model, and (iii) seven Envisat InSAR swaths. Our analysis indicates that the amplitudes of the annual delay vary by up to ~10 cm in this region equivalent to a maximum displacement bias of ~24 cm in InSAR line of sight direction between two epochs (assuming Envisat IS6 beam mode). The stratified tropospheric delay correction mitigates this bias and reduces the scatter due to the stochastic delay. For ~7 years of Envisat acquisitions along the western India plate boundary, the uncertainty of the InSAR velocity field due to the residual stochastic wet delay after stratified tropospheric delay correction using the ERA-I model is in the order of ~2 mm/yr over 100 km and ~4 mm/yr over 400 km. We discuss the implication of the derived uncertainties on the full variance-covariance matrix of the InSAR data

A Network-Based Enhanced Spectral Diversity Approach for TOPS Time-Series Analysis

For multitemporal analysis of synthetic aperture radar (SAR) images acquired with a terrain observation by progressive scan (TOPS) mode, all acquisitions from a given satellite track must be coregistered to a reference coordinate system with accuracies better than 0.001 of a pixel (assuming full SAR resolution) in the azimuth direction. Such a high accuracy can be achieved through geometric coregistration, using precise satellite orbits and a digital elevation model, followed by a refinement step using a time-series analysis of coregistration errors. These errors represent the misregistration between all TOPS acquisitions relative to the reference coordinate system. We develop a workflow to estimate the time series of azimuth misregistration using a network-based enhanced spectral diversity (NESD) approach, in order to reduce the impact of temporal decorrelation on coregistration. Example time series of misregistration inferred for five tracks of Sentinel-1 TOPS acquisitions indicates a maximum relative azimuth misregistration of less than 0.01 of the full azimuth resolution between the TOPS acquisitions in the studied areas. Standard deviation of the estimated misregistration time series for different stacks varies from 1.1e-3 to 2e-3 of the azimuth resolution, equivalent to 1.6-2.8 cm orbital uncertainty in the azimuth direction. These values fall within the 1-sigma orbital uncertainty of the Sentinel-1 orbits and imply that orbital uncertainty is most likely the main source of the constant azimuth misregistration between different TOPS acquisitions. We propagate the uncertainty of individual misregistration estimated with ESD to the misregistration time series estimated with NESD and investigate the different challenges for operationalizing NESD

InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan

We use 2004–2011 Envisat synthetic aperture radar imagery and InSAR time series methods to estimate the contemporary rates of strain accumulation in the Chaman Fault system in Pakistan and Afghanistan. At 29 N we find long-term slip rates of 16 ± 2.3 mm/yr for the Ghazaband Fault and of 8 ± 3.1 mm/yr for the Chaman Fault. This makes the Ghazaband Fault one of the most hazardous faults of the plate boundary zone. We further identify a 340 km long segment displaying aseismic surface creep along the Chaman Fault, with maximum surface creep rate of 8.1 ± 2 mm/yr. The observation that the Chaman Fault accommodates only 30% of the relative plate motion between India and Eurasia implies that the remainder is accommodated south and east of the Katawaz block microplate

InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan

We use 2004–2011 Envisat synthetic aperture radar imagery and InSAR time series methods to estimate the contemporary rates of strain accumulation in the Chaman Fault system in Pakistan and Afghanistan. At 29 N we find long-term slip rates of 16 ± 2.3 mm/yr for the Ghazaband Fault and of 8 ± 3.1 mm/yr for the Chaman Fault. This makes the Ghazaband Fault one of the most hazardous faults of the plate boundary zone. We further identify a 340 km long segment displaying aseismic surface creep along the Chaman Fault, with maximum surface creep rate of 8.1 ± 2 mm/yr. The observation that the Chaman Fault accommodates only 30% of the relative plate motion between India and Eurasia implies that the remainder is accommodated south and east of the Katawaz block microplate

InSAR Time-Series Estimation of the Ionospheric Phase Delay: An Extension of the Split Range-Spectrum Technique

Repeat pass interferometric synthetic aperture radar (InSAR) observations may be significantly impacted by the propagation delay of the microwave signal through the ionosphere, which is commonly referred to as ionospheric delay. The dispersive character of the ionosphere at microwave frequencies allows one to estimate the ionospheric delay from InSAR data through a split range-spectrum technique. Here, we extend the existing split range-spectrum technique to InSAR time-series. We present an algorithm for estimating a time-series of ionospheric phase delay that is useful for correcting InSAR time-series of ground surface displacement or for evaluating the spatial and temporal variations of the ionosphere's total electron content (TEC). Experimental results from stacks of L-band SAR data acquired by the ALOS-1 Japanese satellite show significant ionospheric phase delay equivalent to 2 m of the temporal variation of InSAR time-series along 445 km in Chile, a region at low latitudes where large TEC variations are common. The observed delay is significantly smaller, with a maximum of 10 cm over 160 km, in California. The estimation and correction of ionospheric delay reduces the temporal variation of the InSAR time-series to centimeter levels in Chile. The ionospheric delay correction of the InSAR time-series reveals earthquake-induced ground displacement, which otherwise could not be detected. A comparison with independent GPS time-series demonstrates an order of magnitude reduction in the root mean square difference between GPS and InSAR after correcting for ionospheric delay. The results show that the presented algorithm significantly improves the accuracy of InSAR time-series and should become a routine component of InSAR time-series analysis

Citric acid assisted synthesis of manganese tungstate nanoparticles for simultaneous electrochemical sensing of heavy metal ions

Manganese tungstate nanoparticles (NPs) were successfully synthesized via a simple and cost effective citric acid assisted solution combustion method. Rietveld refinement of PXRD pattern indicates the formation of pure wolframite phase of manganese tungstate. Transmission electron microscope (TEM) analyses confirm the spherical shaped nanoparticles with average particle size ranging from 12 to 26 nm. The ultra-light combustion derived manganese tungstate NPs exhibits mesoporous nature with the estimated band gap 2.9 eV. The synthesized NPs were found to be an excellent electrode modifier for the simultaneous detection of Pb(II), Cd(II), Cu (II) and Hg(II) ions at nanomolar (nM) concentration. The present electrochemical sensor exhibits wide linearity of 10 - 500 nM with limit of detection based on 3 sigma method is 3.3, 3.5, 2.9 and 3.1 nM for Pb(II), Cd(II), Cu(II) and Hg(II) respectively. Additionally, room temperature photoluminescence emission spectrum exhibits an intense band at 445 nm that can be attributed to (1)A(1), ground-state to the high vibration T-1(2) energy levels of distorted WO6] octahedral groups